The surface-state of the topological insulator Bi 2 Se 3 revealed by cyclotron resonance. 
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To date transport measurements of topological insulators have been dominated by the conductivity 
of the bulk, leading to substantial difficulties in resolving the properties of the surface. To this 
end, we use high magnetic field, rf- and microwave-spectroscopy to selectively couple to the surface 
conductivity of Bi2Se3 at high frequency. In the frequency range of a few GHz we observe a crossover 
from quantum oscillations indicative of a small 3D Fermi surface, to cyclotron resonance indicative 
of a 2D surface state. 



PACS numbers: 73.20.-r, 71.18.+y, 71.70.Di, 78.70.Gq. 

A massless Dirac fermion is characterized by a linear 
dispersion, best described in a relativistic framework and 
the discovery of topological insulators which support such 
quasiparticles has sparked substantial interest for both 
fundamental and technological reasons [IH5]. In the last 
year, a large number of surface-sensitive probes have re- 
ported the existence of Dirac quasiparticles, similar to 
those reported in graphene [IHZ] on the surface of sin- 
gle crystals of Bi 2 Se 3 and related compounds [8HTT]. An 
experimental signature of the linear dispersion is that 
the cyclotron resonance (resonant optical absorption in 
the presence of a magnetic field) occurs between Landau 
quantized energy levels E n determined by the expression 

a 

E n = v F ^2(n + j)tieB, (1) 

where v-p is the Fermi velocity, n is the Landau level 
index, 7 is the Berry phase (which depends on the 
topology of the Fermi surface) and B is the magnetic 
field. While recent scanning-probe microscopy experi- 
ments have observed evidence for this unique quantiza- 
tion [T2I [13], cyclotron resonance experiments [14] and 
recent far-infrared measurements |15l 116] on the same 
compounds have only detected signal from the bulk. The 
complications associated with the presence of a bulk 
Fermi surface make transport and cyclotron resonance 
detection of the surface state significantly more difficult 
than in mono-layer graphene |17j . 

In this letter, we report experiments which deconvolve 
bulk and surface properties of Bi 2 Se3, by employing a 
high frequency contactless conductivity techniques. By 
increasing the frequency at which the measurement is 
performed we reduce the skin depth and decrease our 
sensitivity to the bulk conduction electrons. We observe 
2D cyclotron resonance phenomena, indicative of the sur- 
face state coexisting with the bulk 3D Fermi surface, 
as anticipated from our earlier comparison of ARPES 
and quantum oscillation measurements |17) . We find 
that the cyclotron mass of the surface states is enhanced 
over that anticipated from photoemission measurements 



of the Fermi velocity [5], suggesting significant many- 
body renormalization of the surface Dirac fermions of 
this topological insulator. 
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FIG. 1: The magneto-conductivity of Bi2Se3 measured at 
three very different frequencies, top: 4-wire resistance mea- 
sured at 312 Hz, middle: rf-skindepth (frequency shift of res- 
onator) measured at 275 MHz and bottom: microwave cavity 
transmission measured at 71 GHz. All measurements were 
performed at a temperature of 1.5 K with the magnetic field 
applied parallel to the crystallographic c-axis Note that the 
r.f. and microwave measurements were made on the same 
sample, the resistance [17] was measured on a different sam- 
ple from the same growth. The microwave absorption at 2T 
is an isotropic cavity background. The angular dependence of 
the A.C. and microwave measurements are shown in Figure[2] 

Single crystals of Bi2Se3 were grown by slow cool- 
ing a binary melt, conditions are described elsewhere 
|17) . yielding samples with typical dimensions of 2 x 1 x 
0.1 mm 3 . The trigonal c-axis is perpendicular to the 
cleavage plane of the crystals. Low magnetic field Hall 
measurements were used to determine that the bulk car- 
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rier concentrations were 2 ± 0.5 x 10 17 cm" 3 [17 . The 
resistance was measured using standard 4-wire ac trans- 
port at a frequency of 312 Hz. 

The magneto-conductivity was measured using rf- and 
microwave techniques. The former was achieved by in- 
ductively perturbing a resonant tank circuit [TH]. The 
resonant frequency of the circuit is dependent upon the 
inductance of the sample coil, such that as the sam- 
ple skin depth increases with resistance, the resonant 
frequency decreases. For the highest frequency circuit, 
275 MHz, the skin depth is estimated from the bulk con- 
ductivity to be 30 /im, a significant fraction of the sam- 
ple thickness. The conductivity at microwave frequen- 
cies was measured using a cavity perturbation technique, 
whereby the change in cavity transmission at resonance 
reflects the field induced changes in the complex con- 
ductivity of the sample. Two microwave cavities were 
used, one a fixed-angle multi-moded cylindrical cavity 
in the frequency range of 10-40 GHz [TS], the other a 
mono-moded cavity resonating at 71 GHz, that can be 
rotated with respect to the applied magnetic field at cryo- 
genic temperatures [301 |5T]. Both were measured us- 
ing an MVNA spectrometer manufactured by AB-mm. 
For 71 GHz the bulk conductivity, yields a skin depth of 
2 /im, only 50 times the thickness of the depletion region 
estimated for this carrier concentration [17 . It should 
be further noted that these estimates do not account for 
the possible screening effect of a high mobility surface 
state. The oscillating magnetic field in the coil pertur- 
bation technique is parallel to applied magnetic field and 
perpendicular to the (001) crystal surface so only induces 
screening currents in this plane. For both microwave cav- 
ity geometries the sample is located with the oscillating 
magnetic field in the plane of the (001) surface so in- 
duces screening currents both in and perpendicular to 
this plane. Magnetic fields were provided by both su- 
perconducting solenoids (up to 17 T) and by resistive 
Bitter magnets (up to 35 T) at the NHMFL in Tallahas- 
see. Standard 4 He techniques were employed to regulate 
temperature down to 1.5 K. 

Figure[l]illustrates the magneto-conductivity of Bi 2 Se3 
using the three different techniques. The magneto- 
conductivity measured at frequencies below 275 MHz 
is dominated by quantum-oscillations of the bulk Fermi 
surface. The period of the quantum oscillations (in in- 
verse field) is proportional to the Fermi surface cross- 
section, related by the Onsager equation Ak = (n + 
t)^^ [221 [23] . From the angle dependence of the fre- 
quency of the Shubnikov-de Haas oscillations (Figure [2]) 
we deduce that the bulk Fermi surface consists of a sin- 
gle closed ellipsoid, with ellipticity< 2 [T7J [55] . For the 
same sample measured at a frequency of 71 GHz however, 
the signal is dominated by a broad absorption centered 
around 8 T. 

For a different sample from the same growth, Fig- 
ure km) plots the transmission measured at different mi- 
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FIG. 2: The angular dependence of the quantum oscillation 
frequency (measured at 312 Hz) showing slightly anisotropic 
3D behavior (upper points) and the angular dependence of 
the cyclotron resonance field (measured at 71 GHz) show- 
ing 2D behavior (lowerpoints). Measurements correspond to 
the field sweeps shown in Figure [T] Both dashed lines are 
Bo/cos(9), as expected for a 2D phenomenon. The solid line is 
Bo I \Jcos 2 {6) + r] 2 sin 2 (8), corresponding to a 3D Fermi sur- 
face, with an ellipticity, rj — 1.50. 

crowave frequencies, indicating that this is a resonant 
phenomenon with linear frequency field scaling. The 
wide field range over which data from the sample in Fig- 
ure [4] was collected leads us to identify the resonances 
as a single harmonic series, the weaker, lower field res- 
onances occurring at magnetic fields close to ^, | and 
I that of the fundamental. The frequency field scal- 
ing of the fundamental cyclotron resonance condition is 
2.7 GHzT -1 , providing compelling evidence that this is 
not cyclotron resonance originating from quasi-particles 
with an effective mass the same as the conduction band, 
rn* = 0.12 m e , [27] which would lead to a resonance con- 
dition corresponding to 233 GHzT -1 , around 2 orders 
of magnitude greater than our observation. Neither can 
this be a cyclotron resonance harmonic lj = lu c , from a 
band with a mass of 0.12 m e , as at 71 GHz one would 
expect subsequent harmonics to be less than 1 T apart 
in field. This resonance therefore does not originate from 
the bulk band structure, and the remainder of this study 
is dedicated to the investigation of these high-frequency 
features. 

Figures [2] and [4] show the angle dependence of the high- 
frequency resonance for two different samples, where the 
angle 6 is defined as that between the normal to the cleav- 
age plane (the [001] crystal axis) and the field direction. 
The resonance field depends only on the component of 
the field which is perpendicular to the cleavage plane and 
scales as l/cos(9). This is an unambiguous signature of 
a 2D state. As this property only emerges at our highest 
frequencies where the skin depth is smallest, it is natu- 
ral to attribute the origin of this signal to the 2D Dirac 
fermions, known to exist on the surface of the material 
from photoemission experiments [HI HZ] ■ 
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FIG. 3: a) Magnetic field dependence of the cavity transmis- 
sion for a selection of frequencies, top to bottom: 36.3, 27.1, 
18.2, 15.1, 11.4 GHz, offset for clarity, b) The corresponding 
frequency vs magnetic field for the two main cyclotron reso- 
nances observed in the field range 10 to 40 GHz. It should 
also be noted that the two resonances observed do not have 
a simple harmonic relation leading us to believe that they 
originate from opposite sample surfaces. The factor of ss 2 
variation in uc and hence Ef between surfaces is consistent 
with the ARPES results from this growth reported in [T7], see 
surface doping arguments below. 



As we show below the observed frequency scaling, that 
is linear in B rather than yB, as perhaps expected from 
Equation [I] is consistent with the Fermi energy relative 
to the Dirac point being large in comparison to the tem- 
perature and measurement frequency. To model this be- 
havior, and confirm that the linear regime is expected 
for these samples, the absorption, a, is approximated as 
the sum of Lorentzian lines centered upon the frequency 
corresponding to adjacent Landau level separation. The 
intensity for each transition is weighted by the proba- 
bility of occupation of the initial and final states, thus 
ensuring only transitions from occupied to unoccupied 
states contribute. 



a(B) = 



— T - ^— F{E n )[l-F{E n+1 )\ 

(huj - E n - E n+ i) z + t a \ 

(2) 

where F(E) is the Fermi-Dirac distribution function and 
E n is given by equation [T] The magnetic field corre- 
sponding to the fundamental cyclotron resonance condi- 
tion was calculated as the maximum in absorption for a 
given frequency. The calculation was repeated for a range 
of frequencies and for higher harmonics corresponding to 
An 7^ ±1. For absorption line widths, r, less than the 
experimental temperature the result of the simulation 
was found independent of the value used and for a large 
Fermi energy, the resonance condition is independent of 
the value for the additional Berry phase parameter, 7. 
Unlike the case for massive quasi-particles where the gra- 
dient of the resonant frequency/magnetic field condition 
is determined solely by the effective mass, in this case it 




10 20 30 

Magnetic Field (T) 





Angle (°) 



FIG. 4: The angular dependence of the cyclotron resonance 
field and harmonics measured at 71 GHz. Open symbols are 
data points measured in Los Alamos, solid symbols are a sub- 
sequent measurement on the same sample in Tallahassee, note 
the systematic shift to lower field with time. The lines are fits 
to Bo/cos(9) indicating the 2D nature of the resonance. The 
inset shows examples of transmission data offset for clarity, 
top to bottom is 90° to 0° in 10° steps. Note that the reduced 
amplitude at high angles is due to the angular dependence of 
the cavity coupling. 



is determined by both the Fermi-velocity and -energy, 



eBv\ 



(3) 



The simulation yields the same result as evaluating the 
cyclotron mass [23] 



m c (£ , F 



dAk 
2tt dE 



(4) 



at a constant energy (-Ef) as opposed to the non-constant 
energy limit that yields the VB field dependence of the 
cyclotron frequency for band fillings close to the Dirac 
point. This case is satisfied in graphene[5H7], but in the 
present case with high filling, we expect the linear B 
dependence of Equation [3j 

ARPES [HI El [TT] measures the velocity associated with 
the surface state to be 4.2 x 10 5 ms _1 and puts the Fermi 
energy close to the bulk conduction band edge, i.e. ap- 
proximately 200 meV above the Dirac point. However, 
the Fermi surface of the Dirac fermion is set by surface 
charge accumulation which depends sensitively on the 
atmospheric exposure of the sample [IT] . This leads to 
sample and exposure dependence of the 2D cyclotron res- 
onance which depends on the Fermi energy (Equation |3j, 
as shown in Figure [4] and by the different cyclotron res- 
onance conditions observed for the 3 different samples 
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in Figures [T] [3] and [4| In contrast, this behavior would 
not be expected from a bulk cyclotron resonance, which 
depends only on the effective mass. 

Using Ep = 200 meV requires a Fermi velocity of the 
Dirac cone nearly 7 times smaller than that measured 
by ARPES, suggesting substantial many-body renormal- 
ization of the topological surface state Dirac fermions. 
Recent far-infrared reflection measurements |15j indi- 
cate strong electron-phonon coupling, consistent with the 
large dielectric constant [IB]. The authors also report 
coupling between these optic phonons in the 10 meV 
range and magnetic field [TB] . Although electron-electron 
interactions do not usually renormalize the cyclotron 
mass, it has been reported |21j that the thermodynamic 
effective mass, that is probed by the temperature de- 
pendance of quantum oscillations |21j . can be exceeded 
by the cyclotron mass. The observation of an electron- 
electron enhanced cyclotron mass is also in agreement 
with Hubbard-model calculations for interacting systems 
[2~4"j . Recent predictions for the surface state collec- 
tive mode in topological insulators [25J and in particular 
Bi2Se3 indicate that not only does the spin wave couple 
to sound propagation but that Coulomb interactions can 
couple the spin wave and surface plasmon modes. For the 
region of their dispersion where they coexist, it is reason- 
able to speculate whether self-energy corrections due to 
excitation of this spin-plasmon will enhance the observed 
cyclotron mass over that of the bare Dirac dispersion. 
These measurements hence place important constraints 
upon the magnitude of these effects. 

In conclusion, by probing the conductivity at reduced 
skin depths, we have observed a 2D cyclotron resonance 
from a material whose bulk Fermi-surface is 3D. The 
frequency-magnetic held scaling of this resonance is in- 
consistent with the bulk effective mass, but more consis- 
tent with the dispersion and band filling of a Dirac-like 
surface state as observed by ARPES [T7], with substan- 
tial many-body renormalization. 
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